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1.  RESEARCH  SUMMARY 

1.1  Long-Term  Goals 

The  overall  goal  of  this  project  is  to  determine  the  feasibility  of  in  vivo,  non-invasive 
measurement  of  the  complex  elastic  moduli  of  cetacean  head  tissue.  If  this  objective  is  met, 
measurement  systems  could  be  developed  capable  of  measuring  the  complex  elastic  moduli  of 
the  head  tissue  of  live,  stranded  cetaceans. 

1.2  Objectives 

The  technical  objective  is  to  remotely  generate  and  detect  mid-frequency  (1  to  10kHz)  elastic 
waves  within  the  body  of  a  living  cetacean,  using  ultrasound  and  to  use  the  measured 
propagation  parameters  of  these  waves  to  obtain  the  complex  elastic  moduli  by  inversion.  A 
further  technical  objective  is  to  extract  tissue  moduli  in  this  manner  intracranially  This 
objective  carries  considerably  more  technical  risk  since  both  the  wave-generating  ultrasound 
and  the  probe  ultrasound  will  be  attenuated,  distorted  and  scattered  by  the  passage  through 
the  skull. 

1 .3  Approach 

The  approach  is  to  measure  the  complex  shear  and  bulk  modulus,  from  which  all  other  moduli 
can  be  calculated.  The  shear  modulus  is  determined  by  measuring  the  speed  and  attenuation 
of  shear  waves  generated  within  the  tissue  using  focused  ultrasound  as  a  remote  localized 
force  generator.  This  general  approach  to  determining  the  complex  modulus  is  an  application 
of  a  new  medical  imaging  technology  called  elastography.  The  methods  described  by 
Greenleaf  (Chen  et  al,  2002)  and  Sarvazyan  (Sarvazyan  et  al,  1998)  provide  the  basis  for 
shear  wave  generation.  Displacements  are  generated  remotely  in  a  tissue  volume  using  one  or 
more  focused  ultrasound  beams.  In  the  single-beam  case,  the  ultrasonic  carrier  signal  is 
modulated  at  a  low  frequency  Af.  In  the  dual-beam  confocal  configuration,  the  two  ultrasonic 
signal  carrier  frequencies  are  offset  from  each  other  by  2Af.  In  both  cases,  a  radiation  force  at 
the  focal  point  primarily  generates  shear  waves  at  the  frequency  2Af. 
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The  particle  displacements  resulting  from  the  passage  of  the  propagating  shear  wave  are 
detected  remotely  using  a  modified  version  of  an  ultrasonic  Doppler  vibration  measurement 
system  called  NIVAMS  developed  at  Georgia  Tech  (see  Cox  and  Rogers,  1987,  Martin  et  al, 
2002).  The  system  was  modified  to  make  the  transmitting  and  receiving  beams  coaxial  in 
order  to  permit  operation  in  tight  quarters.  Algorithms  have  been  developed  to  enable  the 
magnitude  and  phase  of  vibration  to  be  determined,  as  well  as  the  range  at  which  the  vibration 
is  being  measured.  By  measuring  the  amplitude  and  arrival  time  of  the  shear  wave  at  two 
different  points  the  propagation  speed  and  attenuation  can  be  determined. 

2.  RESULTS 

2.1  Summary 

A  confocal  ultrasonic  elastography  system  was  designed,  developed,  and  tested  with  discrete 
scatterers  and  tissue  phantoms.  The  former  experiments  demonstrated  the  ability  to  use 
ultrasound  for  remote  generation  and  detection  of  vibrations  in  the  audio  frequency  range. 

The  latter  experiments  demonstrated  the  ability  to  remotely  detect  shear  wave  motion 
amplitude  and  phase  as  a  function  of  depth  along  the  ultrasonic  beam  and  determine  the  wave 
propagation  speed.  Results  from  these  tests  were  in  good  agreement  with  results  obtained 
from  a  resident  dynamic  mechanical  testing  facility.  Considerable  effort  was  applied  towards 
refining  the  detection  system  hardware  and  signal  processing  in  order  to  improve  overall 
measurement  quality  and  reduce  noise,  with  the  ultimate  goal  of  reducing  ultrasonic  drive  level 
requirements.  System  simulation  models  were  developed  to  assist  the  vibration  measurement 
algorithm  development  effort  Work  performed  under  this  grant  set  the  stage  for  testing  on 
tursiops  truncatus  tissue  samples  and  live  animals  in  FY09  under  a  follow-on  grant 

2.2  Range  discrimination  in  continuous-wave  ultrasonic  vibrometry 

The  following  is  adapted  from  a  manuscript  in  preparation  for  submission  to  the  Journal  of  the 
Acoustical  Society  of  America. 

2.1.  Introduction 

In  previous  work,  two  fundamentally  different  classes  of  techniques  have  been  used  for 
ultrasonic  vibrometry  One  of  these  involves  a  continuous  ultrasonic  interrogation  signal  and 
the  other  involves  a  pulsed  signal. 
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In  the  first  class  of  techniques,  a  pure-tone  ultrasonic  carrier  signal  is  reflected  from  a 
surface  or  region  of  interest.  A  phase  modulation,  produced  by  the  low-frequency  (LF) 
displacement  of  that  surface  projected  onto  the  carrier  signal’s  propagation  path,  can  be 
extracted  from  the  reflected  signal  and  used  to  compute  the  surface  displacement  In  this 
technique,  which  relies  on  focused  transmit  and  receive  beams  to  achieve  spatial  resolution, 
lateral  resolution  derives  from  the  narrowness  of  the  beams  and  resolution  in  range  derives 
from  the  reflectivity  of  the  target  and/or  the  depth  of  the  intersecting  region  of  the  transmit  and 
receive  beams.  This  technique  was  first  reported  by  Cox  and  Rogers  in  1987'.  In  that  work, 
the  signal  processing  involved  measuring  relative  side-lobe  levels  on  a  spectrum  analyzer  to 
determine  the  amplitude  of  a  pure-tone  LF  vibration  The  technique  was  later  expanded  by 
Martin  et  al"  and  by  Finneran  and  Hastings'"  to  include  LF  phase  measurement  and  real-time 
measurement  of  broad-band  LF  vibration  by  a  variety  of  different  demodulation  techniques. 
These  included  passive  mixers'",  active  analog  multipliers",  phase-locked  loops'v  and  direct 
digitization  of  the  modulated  carrier  signal*.  Pure-tone  ultrasonic  vibrometry  has  been 
implemented  in  both  air-coupled  and  water-coupled  systems  for  a  variety  of  applications 
including  bioacoustic  and  biomedical  imaging,  structural  acoustics,  and  shallow  seismology. Vl 

Pure-tone  carrier  signals  permit  the  transduction  of  small  displacements  (amplitudes  on  the 
order  of  a  1  x  1 0'6  ultrasonic  wavelengths  per  VHz)  but  do  not  allow  for  differentiation  between 
individual  scatterers  contributing  to  the  received  signal  other  than  by  isolating  them  within  the 
region  where  the  transmit  and  receive  beams  overlap,  as  there  is  no  way  to  distinguish  the 
particular  source  of  either  the  modulation  or  the  scattered  energy  for  a  pure  tone.  This 
limitation  also  gives  rise  to  problems  in  sensing  the  vibration  of  surfaces  with  roughness  on  a 
scale  comparable  to  the  ultrasonic  wavelength  of  the  interrogating  signal.  Here,  although  the 
motion  may  be  completely  characterized  by  a  single  time  history,  interference  can  drive  down 
the  energy  in  the  reflected  ultrasonic  signal  and  render  the  extraction  of  this  information 
impossible.v" 

The  pure-tone  technique  is  closely  analogous  to  laser  Doppler  vibrometry,  although  the 
useful  applications  and  signal  processing  constraints  are  somewhat  different.  Although  they  do 
not  penetrate  many  media  of  interest,  optical  signals  can  be  propagated  over  long  distances  in 
air  or  vacuum  without  significant  energy  loss  due  to  spreading.  These  are  therefore  more 
useful  than  ultrasound  when  the  target  is  at  a  great  distance  and  only  its  surface  motion  is  of 
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concern.  Optical  wavelengths  are  typically  much  shorter  than  acoustic  wavelengths  and 
therefore  have  larger  modulation  for  a  given  displacement.  This  simplifies  post-processing 
Unfortunately,  the  much  higher  frequencies  of  optical  signals  constrain  the  available 
techniques  for  extracting  the  modulation  and  may  require  a  mix  down  to  an  intermediate 
frequency  prior  to  demodulation.vl"  Where  the  static  distance  to  the  surface  of  interest  and  its 
reflectivity  are  known  a  priori  to  high  precision,  extraordinarily  small  modulations  have  been 
extracted  from  laser-based  systems  (amplitudes  on  the  order  of  1*10'12  optical  wavelengths), 
which  makes  these  particularly  well  suited  to  characterizing  time  varying  force  fields. IX 

In  the  second  class  of  ultrasonic  vibrometry  techniques,  transmit  signals  comprised  of  a 
series  of  pulses  are  used  to  illuminate  a  volume  of  interest.  The  motion  of  a  particular  target  or 
target  region  within  the  volume  is  then  determined  from  repeated  correlations  of  the 
corresponding  echoed  pulses  with  an  initially  sampled  record  of  either  the  transmitted  or 
received  signal.  These  techniques  permit  different  targets  to  be  simultaneously  interrogated 
and  discriminated  from  each  other  with  a  single  carrier  signal  that  produces  several  time 
separable  echoes.  Unfortunately,  they  require  that  the  motion  of  the  targets  be  relatively  large 
(amplitudes  on  the  order  of  a  tenth  of  an  ultrasonic  wavelength)  because  the  large  crest  factor 
(the  ratio  of  the  signal’s  peak  value  to  its  RMS  value)  and  high  frequency  content  of  the 
transmitted  signals,  along  with  the  nature  of  the  post  processing,  degrade  the  available  signal- 
to-noise.  These  techniques  have  been  developed  primarily  for  biomedical  applications  such  as 
soft  tissue  elastographyx'  XI,X".  The  nearest  non-acoustic  analog  to  these  is  in  pulse-Doppler 
radar  systems*".  These  systems  employ  processing  techniques  which  have  also  been  used  in 
ultrasonic  vibrometryxlv,  although  they  are  primarily  concerned  with  the  measurement  of  nearly 
constant  velocities  rather  than  measuring  vibrations  in  the  audio  frequency  range. 

In  the  present  work,  a  technique  has  been  developed  that  combines  the  sensitivity  of  pure- 
tone  ultrasonic  vibrometry  with  the  range  discrimination  capability  of  pulsed  vibrometry.  The 
new  technique  exploits  a  noise-like  multi-tone  carrier  signal  and  a  novel  post-processing 
algorithm  in  order  to  determine  both  target  range  and  vibration  from  a  continuous  carrier 
signal.  The  technique  is  similar  to  l/Q  (in-phase  /  quadrature)  demodulation.  I/Q  demodulation 
is  generally  applied  to  pure-tone  signals,  but  has  also  been  proposed  to  enhance  correlation- 
based  processors  for  impulsive  signals/'  In  the  present  work,  mixing  is  not  used  to  extract  the 
modulation  as  it  is  traditionally  in  l/Q  demodulation  Instead,  a  modified  pulse-compression 
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method  has  been  developed,  which  offers  several  advantages  over  previously  reported 
techniques.  The  new  technique  was  developed  for  biomedical  use  in  soft-tissue  elastography, 
but  should  be  generally  applicable  to  all  applications  of  ultrasonic  vibrometry. 


2.2.  Basic  Principles 

An  ultrasonic  carrier  signal  is  transmitted  from  and  received  by  effectively  collocated 
transducers  with  narrow  overlapping  beams  directed  along  the  Z  axis.  The  narrowness  of 
these  beams  assures  that  the  problem  is  essentially  one-dimensional  in  Z.  The  carrier  signal  is 
band-limited  and  periodic,  but  is  not  a  pure  tone.  It  is  therefore  made  up  of  a  number  of 
discrete  tones  (A/+1)  distributed  about  a  center  frequency  ( coc )  with  spacing  Aco .  The  carrier 
signal  may  be  exactly  represented  as  the  sum  of  these  components: 


v(')  =  +  mT)  =  Z  v»  cos(<y»/  +  h ) 


(1=0 


m  =  0,1....  M 
N 


[1] 


CO.,  =  CO  - 


-  n  •  A  co 


Where  t  is  time,  T  is  the  period  of  the  signal,  m  is  an  integer,  NAco  is  the  bandwidth,  <f>n\s  the 

phase  required  to  produce  the  waveform,  and  the  total  time  record  is  an  integer  number  of 
signal  cycles  (M)  such  that  0  <  /  <  MT  . 

In  this  section  we  will  assume  the  transmitted  signal  to  be  comprised  of  a  sequence  of 
pulses  that  are  each  discrete  in  time  to  the  extent  that  its  bandwidth  permits  In  the  next 
section  we  will  consider  the  possibility  of  attaining  similar  spatial  resolution  using  pulse 
compression  and  longer  duration  signals  with  identical  bandwidths.  The  pulse-compression 
approach  offers  improved  noise  rejection  at  the  expense  of  increasing  numerical  artifacts  For 
a  specific  scenario,  the  preferred  technique  would  be  determined  by  the  nature  of  exigent 
constraints 

A  single  transmit  cycle  might  be  considered  to  be  a  pulse  with  a  temporal  width  of 

1  2  TC  1 

A /  =  =  —  =  -  T  which  occupies  VN,h  of  a  single-cycle  time  record  and  is  centered  at  a 

liW  NAco  N 
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time  near  the  beginning  of  this  record  given  by/,  =  such  that  the  pulse  is  entirely 

contained  within  the  single-cycle  record  The  received  signal  will  then  be  comprised  of  the 
echoes  resulting  from  the  transmission  of  this  pulse  sampled  between  t=  0  and  t=T.  If  T  is 
selected  to  be  less  than  the  reverberation  time  of  the  medium,  each  pulse  will  sample 

cT 

scattering  from  inhomogeneities  out  to  a  maximum  range  of  zmax  =•—  (where  c  is  the  sound 

speed  in  the  medium),  without  affecting  later  returns  associated  with  the  subsequent  cycle. 
Reverberation  time  is  used  here  and  elsewhere  in  a  generic  sense  as  being  the  time  required 
for  the  returns  to  reach  a  level  of  insignificance  for  the  problem  that  is  being  addressed  rather 
than  a  specifically  defined  parameter  of  the  medium  associated  with  a  level  of  decay.  Thus,  at 
times  greater  than  7",  echoes  should  be  dwarfed  by  the  returns  associated  with  a  subsequent 
cycle.  Since  the  function  V(i)  exists,  there  must  also  exist  a  corresponding  periodic  function 
r(ci)  that  represents  the  transmitted  pressure  in  the  stationary  medium.  Here,  although  it  is  a 
function  of  time,  the  argument  of  P  has  length  dimensions  that  correspond  to  the  one-way 
propagation  of  the  pulse  to  a  point  in  the  medium  Similarly,  the  temporally  periodic  function 
l’(i.ci),  which  accounts  for  the  change  in  pulse  shape  due  to  the  two-way  propagation,  must 
also  exist  as  there  is  no  physical  mechanism  to  disturb  the  periodicity.  The  function  P  also 
includes  absorption  loss  to  and  from  the  scattering  point  z,  as  well  as  diffraction  and  spreading 
effects.  The  received  signal  from  a  stationary  medium  with  scattering  inhomogeneity  may 
therefore  be  represented  as: 

7s(/)=Jj  A(z)P(z,ci  -1z)dz  [2] 

Where  ^(r)  is  the  back-scattering  strength  associated  with  the  point  z.  This  includes  both 
deterministic  (discrete  or  regularly  space  inhomogeneities)  and  random  factors  (speckle 
scattering)  It  accounts  for  the  density,  scattering  strength  of  inhomogeneities  at  z,  the  strength 
of  the  incident  sound  beam,  scattering  directivity  and  the  transduction  coefficients.  It  is  evident 
from  Eq  2  that  the  periodicity  of  P  requires  that  E(i)  is  also  periodic  with  a  period  of  T. 

If  the  medium  within  the  beam  is  vibrating  with  a  displacement  in  the  Z  direction  described 
by  S(=.i)  then  the  received  signal  associated  with  the  m,h  transmitted  cycle  will  have  the  form: 
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/UOC,  =  £  ‘  A(z)P(z,ct  -  2z  -  2S(z,l  -z/c))dz 


[3] 


/>("'  1)7’ 


Note  that  S  is  evaluated  at  t-z/c  since  that  is  the  time  the  sound  reaches  the  scattering  point,  z. 
If  the  displacement  is  required  not  to  have  spectral  components  at  frequencies  other  than 

harmonics  of  ^  and  is  required  to  be  small  with  respect  to  the  acoustic  wavelength  at  ru  ,  then 
rf(c./)will  not  have  the  periodicity  of  P,  and  the  received  signal  can  be  rewritten  as 


M'C;,,,  -  [  ■' A(z)P(z.c,-2z)dz  '  -lf‘A(z)P\z,c, -2ziZS(zJ-zlci^ydz  [4] 

/>0  V - v - ✓ 


modulation : 


ftere  the  evaluation  range  of  the  periodic  functions  has  been  shifted  to  t-0  as  indicated  in  Eq. 
(1)  I”  represents  the  partial  derivative  of  P  with  respect  to  time.  Alternatively,  the  entire  time 
record  can  be  expressed  in  the  form 


t<MT 


n  0 


f>  0 


M)  Z^(-)c°s(^  +  ^,(-)+2 knz\iz  -  2 1  '  YdkJn(z)S{z,()s\n{a}nt  +  i/)l,{z)+2k„z)dz 

-0 


/<A// 


carrier: 


modulation  :  /• 


[5] 


where  ,/t,.(.:)and^„(z)account  for  both  the  original  waveform  shape  and  its  variation  with  range 
in  addition  to  scattering. 

The  two  terms  in  these  expressions  may  be  separated  using  a  comb  filter  because  the  first 
term  is  periodic  and  the  second  is  not.  A  computationally  efficient  way  to  implement  this  filter  is 
to  represent  a  cycle  of  the  carrier  as  the  average  of  multiple  successive  cycles.  In  this  average 
of  the  modulation  should  go  to  zero  as  the  length  of  the  time  record  is  increased  and  the  pass 
bands  of  the  corresponding  comb  filter  are  narrowed.  Thus, 


I<n>r 
£(n>  -l)y 


•iiWK*  |6| 


and,  by  definition 


Emod(t)  =  E(l)-Eiiir(l)  [7] 
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It  is  important  that  the  receive  data  is  sampled  at  an  exact  harmonic  of  the  carrier 
periodicity  to  precisely  preserve  this  periodicity  and  permit  the  separation  of  terms,  as  it  is 
apparent  that  the  modulation  term  is  significantly  smaller  than  the  carrier. 


If  the  width  of  the  pulse  is  sufficiently  narrow  in  time  that  S(z,t )  in  nearly  constant  on  the 
scale  of  c|  <r<c|  /  +  ~  then  the  modulation  may  be  defined  as  a  matrix: 


/::.T  =  -25 
d 


\  2 


+  (m-l)7’]jJ‘/  A{z)P\ct„-2z))dz 


=  m  _  01  A/  ;  n  =  0,1,...— 

2  r 


[8] 


Where  /  is  time  defined  with  regard  to  the  beginning  of  the  cycle  ( 0  <  /„  <T)  and  absolute  time 
is  defined  by t  =  /„  +(m-\)T .  Thus,/,,  is  congruent  to  t  modulo  T,  in  the  same  way  that  time  is 
represented  by  an  analog  clock.  The  matrix ,  although  it  is  two-dimensional,  is  simply  the 

time  vector  Emod(t)  reshaped  into  a  matrix  in  which  the  row  index  designates  time  relative  to 
the  start  of  a  cycle  and  the  column  index  designates  the  relevant  cycle  number  within  the  total 
time  record.  The  new  variable  dn  represents  the  one-way  distance  associated  with  a  round-trip 

delay  of  /„ . 


From  tq.(3)  it  is  apparent  that  the  time  derivative  of  the  carrier  is  given  by 

KoriO  -  c  jl0  A(z)P'(z,ct  -2z)dz  [9] 

Therefore 


5(d,„,r  +{m-\)T)  = 


c) 

'  r*  mod  ^ 
^  n.m 

<  2) 

E'  1 

i « J 

[10] 


where  E'ur\  Motion  of  the  medium  at  pointy,  is  sampled  once  every  T  seconds 

giving  the  motion  at  time  —  ^  —  +  (m- 1  )T.  Thus,  the  received  signals  are  all  essentially 

contained  within  the  first  half  of  each  period.  The  precision  of  the  sampling  time  is  limited  by 
the  pulse  width  At . 
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It  is  apparent  from  Eq.  (4)  that  the  accuracy  of  the  computed  value  of  5{z,t)  will  depend  on 
spatial  variations  in  both  the  true  8(zj)  and  in  the  scattering  term,  a(z).  If  scattering  from  the 
medium  that  is  being  sampled  is  uniform  on  the  scale  of  the  pulse  width,  A/ ,  then  the  spatial 
integration  in  Eq.  (8)  and  Eq.  (9)  can  only  cause  the  computed  displacement  for  a  point,  z,  to 
deviate  from  its  true  value  as  a  result  of  variation  in  S(z,t)  over  A/.  In  contrast,  if  there  is  a 
large  variation  of  scattering  strength  as  a  function  of  range,  A(z),  spatial  misrepresentation  can 
be  caused  by  a  combination  of  the  displacement  and  the  scattering  variations.  This  is  because 
the  relevant  measure  of  pulse  width,  A / ,  depends  on  the  variability  of  the  ultrasonic  scattering. 
A  single  scattering  source  that  produces  an  extremely  large  echo,  can  affect  measurements  at 
significant  distances  because  small  sidelobes  associated  with  its  echo  may  be  larger  than  the 
main  lobes  of  smaller  returns.  Thus  for  practical  purposes,  A/  might  be  considered  to  be  a 
function  of  both  the  bandwidth  and/t(z).  It  may  therefore  also  vary  with  range  in  a  single 

measurement  with  minima  in  regions  of  relatively  uniform  scattering  and  maxima  in  the  vicinity 
of  interfaces.  Because  of  this,  it  is  prudent  to  oversample  the  measured  data  spatially. 

The  measured  displacement  may  be  produced  by  propagating  plane  waves  which  may 
generally  be  described  by: 

8{zj)  =  £ F{z),  cos((Q0  +  /AQ>  +  <p, (z))  [1 1  ] 

U  0 

where  the  highest  significant  frequency  component  of  the  displacement,  Q, ,  is  sufficiently 

small  that  it  can  be  adequately  sampled  at  a  rate  of  ~  samples  per  second  (the  sampling  rate 

inherent  in  the  reshaping  of  £mod(t))  and  lowest  frequency  component  of  interest,  Q(),  may  be 

adequately  sampled  at  1  samples  per  second  (the  lowest  frequency  that  can  be 

MT 

distinguished  by  the  comb  filter).  Although  L  and  AQ  need  not  be  finite,  S(z,t )  will  only  be 

A/  Q. 

represented  in  frequency  bins.  Here  the  characteristic  wavenumber  is  given  by  k,  =  — L, 

2  c, 

where  the  change  in  subscript  for  the  propagation  speed  indicates  that  the  LF  propagation  may 
involve  a  different  mode  than  the  ultrasonic  interrogation  signal  and  that  it  may  involve 
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propagation  in  a  different  direction  (in  which  case  k,  would  be  the  intrinsic  wavenumber 
projected  along  Z,  and  c,  would  be  the  trace  velocity).  Sufficient  spatial  sampling  of  this 


motion  requires  a  pulse  width  of  A/  <  .thus  <  / 

4  (c,k,)  N  4  (c,k,) 

Although  the  spatial  sampling  of  the  measurement  is  dependent  on  the  pulse  width,  the 
spatial  density  of  the  computed  displacement  data,  as  defined  in  Eq.  (10)  is  not,  this  is 
because  the  sampling  rate,  fs ,  of  the  digital  information  must  significantly  exceed  the  pulse 
width  in  order  to  correctly  represent  the  spectra  of  the  carrier  and  modulation  signals  The 

computed  displacement  will  therefore  be  depicted  on  a  spatial  scale  of  A z  =  rather  than 

Ar  L  A/ .  The  more  sparsely  sampled  record  may  be  computed  by  downsampling  the 

computed  displacement  or  the  oversampled  displacement  may  be  interpreted  with 
consideration  of  the  sampling  limits  of  the  problem 


2.3  Pulse  Compression 

In  the  previous  section,  the  transmitted  signal  was  considered  to  be  a  periodic  pulse  This 
is  not  a  requirement  of  the  technique,  and  a  periodic  pulse  may  be  far  from  an  optimal  signal 
for  transmission.  The  reason  for  this  is  that  its  crest  factor  (the  ratio  of  the  signal’s  peak  to  its 
RMS  value)  is  quite  high.  It  therefore  makes  poor  use  of  the  dynamic  range  available  in  the 
components  that  must  comprise  an  actual  vibrometer  system.  The  bandwidth  and  periodicity  of 
the  transmitted  signal  are  its  only  truly  essential  characteristics.  A  pulse  shape  and  width  can 
be  synthesized  from  the  acquired  data  to  suit  the  available  bandwidth  and  the  spatial  sampling 
needs  of  the  problem.  This  is  not  dependent  on  any  temporal  features  of  the  transmit  signal 
other  than  its  bandwidth  and  periodicity. 


Generally,  a  pulse  compression  of  a  signal  C(t)  that  is  the  output  of  a  linear  time-invariant 
system  and  was  originally  generated  by  a  drive  V(t)  using  a  compression-pulse  shape  Vi(t), 
can  be  expressed  as: 


'lDFT{C{t))-DFT{V,(i))Y 

DFT(v(i  j)  J 


=  IDFT 


vM 


[12] 
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where  DFT  and  IDFT  represent  a  discrete  Fourier  transform  and  its  inverse.  This  general 
relationship  is  valid  for  all  pulse  shapes,  Vi(t),  provided  that  they  do  not  span  a  greater 
frequency  band  than  the  drive,  V{t),  as  that  would  give  rise  to  singularities.xvlxv"  Pulse 
compression  via  cross  correlation  (inter  correlation)  of  the  transmitted  and  received  signals, 
which  is  optimal  for  certain  applications  such  as  single-target  detection  against  white  noise,  is 
a  subclass  of  the  compression  described  by  Eq.  (12)  where  the  compression-pulse  shape, 
V,(t),  is  chosen  to  be  the  autocorrelation  of  the  original  drive  signal,  V(t). 

The  data  C(f)  will  have  the  same  spectral  components  as  the  drive  V(t).  In  this  case, 
relationship  between  the  acquired  data,  C(t),  and  its  pulse-compressed  representation,  /(f),  is 
given  by: 

c(0=ZAcos(r<V+<O  =>  l{t)=Y,A»B»cos((0«l  +4>n  +  <P„)  n3i 

11=0  «= 0 


where  B 


|v,H 

|v(4 


and  (pn  =  tan 


imag 


V,  (,„) 

v(<„) 


\\ 


real 


v,M 

v(«) 


For  the  ultrasonic  vibrometer,  a  compression  of  the  carrier,  as  expressed  in  Eq.  (5),  with  a 
pulse  {Vi)  will  add  similar  leading  real  coefficients  and  phase  terms  to  the  expression  yielding: 


1  ur{')=  J  ,  X  An (ZR  cos(co„/  +  <j>„ (z)  +  <pn  +2 k nz)dz  [14] 


n  0 


where  Bn  is  normalized  by  the  original  drive  amplitude  (|v(<yn)j ). 


A  zero  order-hold  (ZOH)  interpolation  of  the  discrete-tone  spectra  of  the  drive  signal,  V(t), 
and  the  compression  signal,  Vi(t),  may  be  defined  such  that:  V',  (ro)=  V, (ry„),  and 

V'(c;j)  V (<u, , )  for  -  ~°  <  (o  <  Mn  +  .  These  spectra  are  a  step-like  construction  for 

computation  that  correspond  to  the  spectral  bin  size  of  records  with  length  T  rather  than  MT 
Another  way  of  viewing  this  construction  is  that  the  cycle-by-cycle  analysis  of  the  received 
signal  described  in  the  previous  section  would  not  distinguish  between  carrier  tones  and  their 
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associated  modulations  (provided  that  those  modulations  are  sufficiently  sampled  at  the 
fundamental  frequency  of  the  carrier),  as  they  will  be  represented  in  a  single  spectral  bin. 

A  compression  of  the  modulation  using  the  ZOH-interpolated  spectra  of  the  /-pulse  will 
yield. 

Anod (')  =  2 f  (*)*.  s\n{eo„i  +  <j)n {z)+2k„z  +  <p„)c/z  [15] 

//=0 

This  compression  is  equivalent  to  the  full  time  record  that  is  recovered  by  compressing  the 
modulation  record  in  sequential  blocks  of  size  T.  It  is  conceptually  useful  to  visualize  the 
process  in  this  way  because  each  7- sized  block  of  the  modulation  record  is  causally  related  to 
a  T-sized  block  of  the  carrier  (just  as  it  was  for  pulsed  signals),  because  7  is  longer  than  the 
reverberation  time  of  the  medium  It  is  also  computationally  useful  to  implement  the 
compression  in  this  way  because  the  received  signal  need  not  be  transformed  into  the 
frequency  domain  in  its  entirety  in  order  for  the  pulse  compression  to  be  performed.  All  of  the 
transforms  described  in  Eq.  (12)  may  be  performed  piece-wise  on  the  modulation,  and  the 
carrier  can  be  distinguished  from  the  modulation  in  the  time  domain  using  Eqs.  (6)  and  (7) 
rather  than  in  the  frequency  domain 

Comparison  of  Eqs. (13)  and  (14)  with  the  corresponding  terms  in  Eqs. (4)  and  (5)  implies 
that  the  compressed  record  of  the  modulation  and  carrier  signals  differ  only  by  the 
compression  coefficients  and  thus  can  be  used  to  discriminate  spatially  to  the  same  extent  as 
if  the  pulse  sequence  were  transmitted  rather  than  synthesized.  This  implication  is  not  entirely 
accurate  in  that  the  sampling  times  for  the  displacement  within  a  cycle  of  the  carrier  signal  are 
dependant  on  both  the  compressed  pulse  shape  and  the  actual  transmit  signal  and  this  can 
lead  to  numerical  artifacts  that  would  not  be  present  if  the  transmit  signal  were  an  actual 
periodic-pulse  signal.  These  artifacts  are  not  intrinsic  to  the  manipulations  associated  with  the 
pulse  compression.  In  section  IV  these  operations  are  shown  to  be  a  convenient  and  artifact- 
free  post  processing  technique  for  data  acquired  using  periodic-pulse  carrier  signals.  The 
benefits  of  a  particular  carrier  signal  over  another  are  related  to  noise.  The  noise  in  the 
compressed  record  is  determined  by  both  the  measurement  noise  and  to  the  transmitted 
signal  Thus,  increasing  the  energy  in  this  signal  will  reduce  the  non-displacement  noise  in  the 
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displacement  computed  from  the  pulse  compressed  record.  This  will  be  further  addressed  in 
section  IV. 

The  principle  issue  for  vibrometry  is  resolving  clutter  and  thereby  correctly  representing  the 
loci  of  vibrations  rather  than  detecting  a  target  against  noise.  A  pulse's  side-lobe  levels  can  be, 
and  often  are,  more  important  than  the  main-lobe  width  in  this  application.  A  procedure  for 
determining  pulse-shape  optimality  is  not  readily  apparent  given  the  various  and  unusual 
constraints  of  this  problem.  As  previously  mentioned,  the  merit  of  a  particular  compression 
pulse  shape,  F,  (r),  depends  on  the  expected  spatial  scale  and  magnitude  of  variation  in 
scattering.  In  experiments,  the  received  signal  measured  at  short  ranges  has  been  dominated 
by  transmitter-receiver  crosstalk,  which  masquerades  as  an  additional  stationary  scattering 
source,  and  by  scattering  from  the  interface  between  the  medium  of  interest  and  the 
transducer  coupling  fluid  These  have  been  the  most  significant  factors  considered  in  selecting 
the  pulse  shape  that  was  used  to  process  the  experimental  data. 

For  the  laboratory  experiment,  the  chosen  compression  pulse  was  a  Gaussian-windowed 

2 

sine  wave  computed  from  V,  =  sin(ry  (/  -/t  }  with«  =  — ^  .  The  pulse  is  depicted  in  Fig. 

1  Although  this  is  functional,  it  is  not  in  any  obvious  way  optimal  for  the  application.  Because 
of  the  band-limited  processing,  the  effective  pulse  shape  is  more  accurately  described  by  this 
pulse  passed  through  an  ideal  boxcar  filter  applied  in  the  frequency  domain.  Because  of  this, 
the  pulse  compression  served  both  to  temporally  compress  the  waveform  and  to  band-pass 
filter  the  already  digitized  data  The  filtered  compression  pulse  includes  some  low-level  ringing 
(both  coda  and  onset)  not  characteristic  of  its  archetypical  form 

The  selected  compression  pulse  is  slightly  wider  than  a  cardinal-sine  pulse  in  main-lobe 
width,  but  has  smaller  sidelobes  at  all  levels  below  about  -10  dB.  This  is  advantageous  given 
the  disparity  between  the  large  crosstalk  and  reflection  from  the  surface  of  the  media  of 
interest  and  the  much  smaller  scattering  from  their  interiors  that  was  observed  laboratory 
experiments.  Although  the  cardinal-sine  pulse  might  seem  optimal  as  it  represents  the 
autocorrelation  of  a  spectrally  flat  transmit  signal  and  therefore  an  optimal  signal  for  single¬ 
target  detection  against  noise,  this  pulse  is  problematic  in  this  application  because  of  the 
number  and  extent  of  its  sidelobes. 
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Another  important  aspect  of  the  compression  pulse  is  that  both  the  pulse  and  its  Hilbert 
transform  are  similarly  discrete  in  time.  The  reason  for  this  is  apparent  in  the  following  section. 


2.4  The  Effects  of  Non-Displacement  Noise 

From  a  transduction  standpoint,  noise  is  only  problematic  if  it  is  not  representative  of 
motion  in  the  region  of  interest.  Motional  noise,  which  may  be  a  problem  in  many  vibrometer 
applications,  should  be  detected  and  properly  represented  by  any  vibrometer  and  may  be 
considered  a  component  of  the  signal  that  is  of  interest.  In  contrast,  non-displacement  noise, 
limits  the  ability  of  the  vibrometer  to  correctly  represent  this  signal.  A  vibrometer’s  insensitivity 
to  this  form  of  noise  is  a  figure  of  its  merit. 

Although  Eq.  (10)  in  conjunction  with  an  appropriate  pulse  compression  would  appear  to 
be  sufficient  for  computing  £(::,/),  it  is  inordinately  susceptible  to  noise.  This  is  because  the 

denominator  (sampled  at  fs)  has  zero  crossings  roughly  every  -  seconds.  At  or  near  these 

points  the  expression  can  blow  up  and  may  either  have  the  form  of  the  displacement  signal 
(because  the  slight  frequency  offset  associated  with  Q,  causes  the  modulation  amplitude  and 

the  carrier  to  pass  through  zero  at  different  points  in  time)  or  it  will  have  the  form  of  noise 
(because  measured  noise  will  be  incoherent  in  the  pulse  compression  and  is  unlikely  to  pass 
through  zero  coincidently  with  the  carrier).  If  a  second  pulse  compression  is  performed  on  the 
same  set  of  data  using  the  Hilbert  transform  of  the  first  compression  pulse  shape  (hereafter 
designated  by  Q  in  contrast  to  /)  then  the  two  compressed  carriers'  zero  crossings  should  not 
be  coincident  unless  there  is  a  null  in  the  scattered  energy  (i.e.  a  true  cancellation  of  the 
signals  from  different  scattering  sources  rather  than  just  a  zero  crossing  in  a  time  waveform) 
This  is  because  the  two  pulses  are  in  quadrature  by  definition  It  is  therefore  advantageous  to 
combine  the  results  of  these  two  pulse-compressed  computations  of  displacement  rather  than 
working  with  one.  Simply  averaging  the  displacement  expressions  together,  would  clearly  be 
inadequate  for  this,  because  the  averages  would  be  dominated  by  the  aforementioned 
singularities.  Instead  of  this,  a  weighted  average  of  the  computed  displacements  was 
considered  using  the  form: 
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[16] 


S(z,t)  = 


rtz)8iM+&)8QM 

${z)+Z(z) 


or 


S(z,t)  = 


f  ( dl  \  (dO  ^  ^ 
A  Z1SSL  +%  U^ar 


V 


\  dt  J 


V  dt  ) 


[17] 


Where  j d  - 


-£ - and  #  = - £- 


(dlcar) 

dQcar_ ' 

l  dt  j 

z 

y  dt  j 

Eq.  (17)  is  generally  valid  for  all  values  of  A  and  ZB  (the  weighting  coefficients)  and  is 

clearly  superior  to  nominally  equivalent  expressions  in  which  the  denominator  is  represented 
as  a  product  rather  than  a  sum  of  carrier  values  (i.e.  carrier-independent  weightings),  as  those 
would  permit  singularities  similar  to  Eq.  (10)  when  either  carrier  term  is  zero  (i.e.  more 
frequently).  With  finite  weighting  coefficients,  Eq.  (17)  only  permits  these  singularities  when 
both  carrier  terms  are  either  zero  or  are  of  equal  magnitude  and  opposite  sign. 

If  the  noise  in  the  received  signal  (Eq.  (3))  is  additive,  then  the  carrier  portion  of  this  signal 
will  have  additive  noise  equivalent  to  the  measured  noise  comb  filtered  at  the  N+*\  carrier 

N 

frequencies  (i.e.  noise n ).  Similarly,  the  modulation  will  have  additive  noise  equivalent  to  the 

rt-0 

N  L 

measured  noise  comb  filtered  at  the  L-(N  + 1)  modulation  frequencies  (i.e.  ^  V  noise n, ). 

n=0  /= 0 

Since  L>  1  and  each  noise  component  is  incoherent,  modulation  noise  would  be  about  Jl 
larger  than  the  carrier  noise  if  the  measurement  noise  were  white.  The  pulse-compressed 
terms  will  also  have  additive  noise,  which  will  be  the  measured  noise  passed  through  the 
appropriate  comb  filter  and  through  a  synthetic  transfer  function  associated  with  the  pulse 

N  ML 

compression.  Thus,  Jfcar  =  ZB,  n°isen  ■  *  ^  mod  =  2X1  noisenl  . 

n=0  n- 0  1=0 

This  noise,  unlike  the  signal,  has  no  spatial  dependence  because  it  does  not  originate  from 
scattering  loci.  It  will,  however  have  an  apparent  spatial  dependence  as  it  is  processed  along 
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with  the  signal.  The  noise  will  alias  from  its  true  ultrasonic  frequencies  into  the  band  of  the 
displacement  records  when  it  is  discretely  sampled  in  intervals  of  T.  Consecutive  points  in 

space  will  represent  consecutive  samplings  of  the  noise  in  intervals  of  —to  create  its  apparent 

f s 

2  dependence.  Since  the  phase  of  the  measured  noise  is  random,  the  phase  change 
associated  with  the  pulse  compression  may  be  neglected.  With  this  noise,  carriers  and 
modulations  will  be: 

/„,(<)=  f'  Y.AXz)B,coicoj+^+<p,+2k,z)iz  +  M-aXl)  (18a) 

n= 0 

Qcar  (0  =  f CT  X  A*  ( Z)Dn  * H(0n‘  +  </>n  +  <Pn  +  2k„z)dz  +  (/))  [1  8b] 

n= 0 

7mod  (0  =  2  f7-  S(Z> 0Z  k  nAn  C Z)Bn  S'nK'  +  4,  (*)  +  2KZ  +  <P  n  ¥Z  +  Anod  (0  t1  8C] 

n- 0 

^mod (0  =  2 1  T <*(z>0ZM, (ZR cos(ay  +  (j)n (z)+  2 knz  +  <pn)dz  +  (/))  [1 8d] 

n= 0 


where  K  denotes  a  Hilbert  transform.  It  is  necessary  to  preserve  the  phase  associated  with 


this  transform  because  the  two  noise  components  in  Eqs.  (18a)  and  (18c)  are  in  quadrature 
with  the  noise  terms  in  Eqs.  (18b)  and  (18d)  by  definition.  The  displacement  has  been 


assumed  small  in  the  analysis  and  therefore  ,  Qm0i(t)<<  — 


dt 


dt 


thus,  with  the 


additive  noise,  Eq.  (17)  has  the  form 

/«-  iZAA  + 


#(ZA 


apparent 


d AdQ- 


[19] 


dt 


/ 


dt 


'zj 


S(zj) 


It  is  clear  from  this  equation  that  an  appropriate  selection  of  A  and  S  will  improve  the  noise 

immunity  of  the  measurement  since  both  the  compressed-modulation  noise  and  its  Hilbert 
transform  have  equal  energy,  but  the  two  components  of  the  modulation  may  be  widely 
disparate  in  energy  content  depending  on  the  values  of  the  compressed  carrier  signals 
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associated  with  the  location,  z.  If  this  expression  is  considered  for  a  point  in  the  time  record 
where  /mrf(z,r)=£Gmd(z,<)  and 


dl  ^ 

car 

—  c 

dt  ) 

—  G 

dt  ) 

» 

?r.,1  pn, 

°  \Z’  1  /apparent  =  0[Z,t )  + - - - y~  ^ -  l^OJ 

2  dQ™ 

Ch  \  dt 


The  modified  noise  term  NmoXzd)  has  both  an  apparent  spatial  dependence  and  an 

apparent  low-frequency  time  dependence  due  to  downsampling  and  aliasing  in  the  reshaping 
of  the  modulation  vectors  into  matrices.  The  resulting  apparent  displacement  noise  preserves 
the  quadrature  Hilbert-transform  relationship  because  the  aliased  terms  represent  identical 
mixes  (i.e.  downshifts)  for  both  components  of  the  numerator  in  the  displacement-noise  term.  If 

A 

we  define  X  = — (the  ratio  of  the  shadings),  then  noise  energy  (represented  here  in  the  low- 
frequency  domain)  is  minimized  by  minimizing  the  expression: 


X  +  i 
Xs  + 1 


[21] 


therefore 

X2£  +  X(\-£2)-e  =  0  [22] 
The  roots  of  this  expression  are  at: 


Xx=£ 

x  =_I  P3] 

£ 


Where  the  1st  root  represents  the  noise  minimum  and  the  second  is  a  singularity.  This 


defines  the  optimal  weighting  as  A  =  A(z)  = 


dL 


dt 


and  =  rB{z)  = 


dQa 


dt 


Thus,  the  optimal 


weighting  coefficients  in  Eq.  (16)  (g  and  £)  are  proportional  to  instantaneous  energy  in  the 
compressed  carrier  signals.  Substituting  these  values  yields: 
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[24] 


S(z,t)  = 


*Q^Md^L 


( , 


dL 


dt 


+ 


J  2 


dQcar 

dt 


>2  \ 


If  noise  is  included  in  the  exact  expression,  then 


d{z,t)  =  S{z,t)+ 


dt 


dt 


dI  cm 

dt 


)z 


dQ„ 


dt 


[25] 


At  each  point  in  the  apparent  displacement  time  record,  the  downsampled  modulation  noise 
has  nominally  the  same  power,  however  the  apparent  displacement  noise  is  diminished  by  the 
energy  in  the  reflected  carrier.  Thus  the  apparent  displacement  noise  in  Eq.  (25)  will  be  high 
for  regions  of  apparently  weak  scattering  and  low  for  regions  where  the  scattering  perceived  by 
the  receiving  transducer  (i.e.  accounting  for  attenuation,  shadowing  etc.)  is  strong.  Non¬ 
displacement  noise  should  not  otherwise  vary  with  range  as  the  ultrasonic  noise,  unlike 
vibrational  displacement,  is  not  associated  with  any  particular  scattering  loci.  The  shading 
represented  in  Eq.  (24)  may  seem  obvious  given  the  aforementioned  need  to  avoid  nulls  in  the 
denominator,  however  it  should  be  noted  that  the  optimality  of  this  form  applies  to  all  values  of 
s  whereas  nulls  due  to  zero-crossings  in  Eq.  (10)  exist  by  definition  only  where  t  e(0,oo). 
Even  for  values  of  z  at  which  one  of  the  compressed  carriers  is  significantly  larger  than  the 
other,  there  is  still  an  advantage,  albeit  small,  to  including  the  result  of  both  compressions  in 
the  computation.  Where  both  compressed  carriers  are  equal  in  amplitude,  this  advantage  is  3 
dB.  It  should  also  be  noted  that  Eq.  (24)  is  not  intrinsically  free  of  singularities  as  it  is  possible, 
but  less  probable,  for  the  denominator  to  be  zero  if  there  is  destructive  interference  in  the 
scattering  from  multiple  targets  contributing  to  the  signal  associated  with  a  particular  range,  z 
This  is  a  different  scenario  than  simply  crossing  zero  in  a  time  waveform.  Here,  it  is 
advantageous  to  downsample  spatially  by  representing  the  displacement  at  a  point  as  the 
weighted  average  of  displacements  from  several  adjoining  points,  where  the  weighting 
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coefficient  for  each  term  is  the  corresponding  denominator  from  Eq.  (24).  This  downsampling 
is  not  intrinsically  synonymous  with  downsampling  to  account  for  the  difference  between  the 
compressed  pulse  width  and  the  sampling  rate,  and  thus  does  not  obviate  the  need  to  consider 
the  compression-pulse  width  in  analyzing  the  data  regardless  of  the  existence  or  absence  of 
nulls 

In  both  numerical  simulations  and  in  experiments,  Eq.  (24)  can  be  shown  to  offer  superior 
noise  immunity  to  a  single  pulse  compression.  This  is  shown  in  Fig  2  where  the  ratio  of  the 
apparent  displacement  energy  in  signal-free  noise  processed  using  Eq.  (10)  to  the  same 
record  processed  using  Eq  (24)  is  depicted  as  a  histogram.  Here  two  cases  are  shown.  In  the 
first  case,  (A),  the  record  is  a  synthetic  record  generated  from  two  point  scatterers  that  were 
stationary  (i.e.  signal-free)  with  additive  white  noise.  In  the  second  case,  (B),  the  record  was 
from  a  laboratory  measurement  on  a  tissue  phantom  where  the  scattering  was  effectively 
continuous  and  the  noise  was  from  a  variety  of  sources  some  of  which  are  discussed  in 
subsequent  sections,  but  almost  certainly  combined  both  displacement  noise  and  non- 
displacement  noise  It  is  clear  from  Fig.  2  that  the  use  of  Eq.  (24)  in  lieu  of  Eq.  (10)  reduces 
apparent  displacement  noise  at  all  ranges  although  the  magnitude  of  this  improvement  varies 
considerably. 

Eq.  (24)  is  similar  to  expressions  for  l-Q  demodulation  of  a  pure-tone  ultrasound  signaT 
where  the  range  computed  by  time  of  flight,  z,  replaces  the  physical  location  of  the 
interrogating  ultrasound  transducers  and  relative  phase  is  the  result  of  pulse  compressions 
rather  than  analog  mixes.  Thus,  the  new  signal-processing  technique  is  similar  to  the  point- 
wise  measurement  of  displacements  using  a  pure-tone  carrier  signal  as  transducers  are 
moved  across  a  reflecting  surface.  In  the  case  where  the  carrier  signal  is  a  periodic  pulse, 
referring  to  the  numerical  manipulations  as  a  “compression”  is  a  misnomer.  Here,  the  results  of 
the  two  compressions  need  not  be  more  temporally  discrete  than  the  unprocessed  data.  The 
compressions  may  serve  only  to  shift  and/or  distort  the  waveforms  and  to  generate  a  second 
record  that  is  in  quadrature  with  the  first  in  order  to  reduce  the  effects  of  noise. 

2.5  Method  of  Implementation 

In  an  experiment,  a  transmit  signal  made  up  as  a  sum  of  100  to  400  discrete  tones  with 
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analytic,  random,  or  numerically  optimized  phases  is  used.  A  complete  time  record  of  both  the 
transmitted  and  received  signals  is  digitized  and  saved  with  a  high-speed  digitizer 
synchronized  to  the  signal  source.  The  carrier  is  separated  from  the  modulation  in  the  received 
signal  by  time  averaging  and  subtraction  as  shown  in  Eqs.  (6)  and  (7).  The  carrier  is 
compressed  twice  (/  and  Q  compressions)  as  shown  in  Eq.  (12)  using  previously  synthesized 
compression  pulse  spectra  and  the  recorded  transmit  signal  record.  The  modulation  signal  is 
similarly  compressed  using  the  identical  coefficients  in  successive  blocks  of  one  cycle  length. 
Each  of  these  blocks  is  then  used  to  compute  the  displacements  associated  with  a  particular 
time  (range)  within  that  cycle  using  equation  24.  These  are  stored  in  a  single  column  of  a 
matrix  The  other  columns  of  this  matrix  are  populated  by  the  processing  of  subsequent  blocks 
of  the  modulation  record  Once  populated,  this  matrix  represents  displacements  as  functions  of 
oversampled  range  (columns)  and  time  (rows)  Where  it  is  appropriate,  the  displacement  is 
then  spatially  downsampled  by  a  weighted  average  of  adjoining  points.  Here  a  twenty-point 
average  has  been  typical. 


2.6  Comparison  with  Other  Techniques 

It  is  possible  to  show  that  the  forgoing  amounts  to  a  time-efficient,  noise-immune  pulse 
correlation  technique  ZOH  pulse  compression  of  the  entire  (un-separated)  received  signal 
(without  small  amplitude  assumptions)  will  be  of  the  form. 

/,„  =  [ '  X  A"  (Z)B»  C0S(*V  +  <!>n  +  <P„  +  2 kn  [z  +  £(r,/))}£  [26] 

n~  0 

Cross  correlation  of  the  transmitted  and  received  signals  (Eqs.  (2)  and  (3))  is  the  special 
case  where  Bn  =V  and  <pn  =-</>„.  Here, 


.v  =  V  l  '  Y.An(=)C0S 

/?» o  \ 


2z  2  S{zj)X) 
COJ  t  +  -  + - - - 

c  c  J 


dz  [27] 


For  a  given  LF  time,  the  term  ‘  J  is  a  delay  associated  with  the  return  from  scatterers 

c 

in  the  vicinity  of  z  that  corresponds  to  a  particular  phase  point  on  X rcc .  Usually  only  a  signal 
peak  would  be  considered,  but  the  relationship  is  valid  for  any  defined  signal  phase.  Peaks,  in 
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this  case,  correspond  to  phase  points  at  which  a  quadrature-pulse  compression  is  not 
beneficial  as  its  result  would  be  weighted  by  zero  in  the  evaluation  of  Eq.  (24).  These  points 
are  thus  analogous  to  those  represented  by  the  left-most  bars  of  Fig.  2.  By  tracking  each  of 
these  points  successively  as  a  function  of  t,  it  is  then  possible  to  reconstruct  S(z,l,)  without 

relying  on  the  assumption  of  small-amplitude  displacements  to  simplify  Eq.  (27).  This 
reconstruction  is  nominally  independent  of  the  temporal  distribution  of  energy  in  the  actual 
transmitted  signal,  but  is  dependant  on  the  signal-to-noise  ratio  (SNR)  of  the  received  signal 
Thus  a  transmit  signal  with  a  low  crest  factor  (a  continuous  or  circular  multisine  signal)  should 
provide  a  better  estimate  of  S(z,i)  than  a  sequence  of  impulsive  signals  separated  by  silence, 
if  all  other  factors  were  equal,  because  the  Cramer-Rao  lower  bound  (CRLB)  of  the  estimation 
would  be  reduced  by  the  ratio  of  transmit-signal  energies  in  each  of  the  two  cases. xvmxlx  In  a 
general  sense,  temporal  concentration  of  energy  in  the  actual  transmit  signal  is  neither 
required  nor  generally  desirable  for  vibrometry  similar  to  the  way  that  it  is  not  advantageous  for 
static  echo-location  or  other  linear  system-identification  tasks. xx  In  a  peak-amplitude-limited 
system  with  large  ultrasonic  and  vibrational  bandwidths,  the  improvement  in  sensitivity 
associated  with  a  continuous  rather  than  impulsive  carrier  signal  could  be  substantial  The 
extent  to  which  the  explicit  displacement  expression  in  Eq  (24)  is  more  noise-immune  in 
practice  than  delay  tracking  in  the  compressed  signal  defined  by  Eq.  (27)  is  largely  dependant 
on  the  algorithm  that  is  implemented  to  track  the  delay. 

From  the  standpoint  of  processing  speed,  it  is  clear  that  a  compression  similar  to  the  one 
represented  in  Eq.  (27)  could  be  processed  in  a  manner  similar  to  that  represented  by  Eq  (10) 
and  that  the  susceptibility  of  this  expression  to  noise  could  be  dealt  with  by  only  considering 
those  points  that  correspond  to  maxima  of  the  denominator  as  these  are  equated  with 
correlation  peaks  This  should  yield  a  measurement  sensitivity  that  is  at  least  on  a  par  with  that 
achieved  through  delay  tracking  of  correlation  peaks  as  the  expression  for  displacement  is 
exact  and  the  implicit  zero  weighting  of  the  un-computed  terms  is  optimal.  Computing  the 
displacement  in  this  way  should  require  far  fewer  integral  transforms  than  repeated 
correlations,  and  does  not  require  that  the  measured  data  be  resampled.  The  current 
implementation  of  the  algorithm  has  a  total  run  time  of  a  few  seconds  using  Matlab  7.0  on  an 
older  desktop  PC  to  compute  3000  displacement  time  histories  from  0.5  seconds  of  data 
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sampled  at  6  MHz  (i.e.  3xl06  computed  displacements  in  time  and  range).  Much  of  this  run 
time  is  expended  accessing  and  allocating  memory.  This  speed  is  achieved  despite  the  fact 
that  the  record  lengths  and  signal  frequencies  were  not  chosen  to  make  optimal  use  of  the  FFT 
and  none  of  the  operations  are  performed  using  dedicated  hardware  Pseudo  real-time 
implementation  (processing  delays  equal  to  or  less  than  LF  record  lengths)  should  be  easy  to 
achieve  with  this  technique.  At  the  heart  of  the  algorithm  are  three  FFTs  (one  with  real  input 
and  two  with  complex  input).  If  the  carrier  is  sampled  at  6  MHz  and  processed  in  50-ps 
intervals  (2  KHz  sampling  of  the  displacement)  to  produce  a  two-dimensional  image  from  the 
incoming  data  stream,  as  in  the  current  off-line  implementation,  these  transforms  optimally 
implemented  would  collectively  require  approximately  1  GFLOPS  i.e. 
2 -3  1000 -5 -3000 -log, (3000).  Processing  a  similar  spatial  /  temporal  density  of  data  using 

repeated  correlations  would  require  the  equivalent  of  12xl06FFTs  per  second,  each  of  which 
would  involve  data  resampled  with  a  step  size  comparable  to  the  intended  measurement 
resolution  (i.e  conceptually,  as  high  as  1  THz  resampling  for  the  measurement  of  a  1  nm 
displacement).  Comparative  computational  savings  should  be  realized  when  the  length  of  the 
individual  correlation  records  exceeds  10  resampled  points.  If  it  were  ultimately  possible  to 
achieve  comparable  resolution  with  a  repeated-correlation  algorithm,  this  record  length  would 
likely  be  on  the  order  of  thousands  of  points,  and  the  computational  time  savings  associated 
with  the  alternative  algorithm  should  be  substantial. 

These  speed  and  resolution  improvements  cannot  be  easily  realized  with  a  software 
modification  to  an  existing  pulse-echo  system  because,  even  with  a  periodic  impulsive  transmit 
signal,  implementation  of  the  algorithm  requires  that  the  trigger  timing  of  the  transmit  pulse  be 
extraordinarily  stable  (on  the  order  of  the  clock  jitter  of  the  digitizer)  and  the  pulse  shape  be 
extremely  repeatable  (on  the  order  of  one  part  in  105).  Here  the  analogous  correlation-based 
system  would  require  correlating  a  sequence  of  receptions  from  a  particular  scatterer  with  a 
single  record  of  the  transmit  pulse.  Any  jitter  in  transmit-pulse  timing  or  shape  will  then 
masquerade  as  displacement  noise  If  the  processing  algorithm  would  otherwise  be  capable  of 
achieving  the  CRLB  associated  with  the  data  record,  it  would  be  degraded  by  this  false 
amplification  of  the  noise  floor  and  thus  the  CRLB.  An  algorithm  that  exploits  the  transmit-pulse 
record  specifically  associated  with  a  received  signal  in  order  to  avoid  the  jitter  in  trigger  timing 
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would  also  be  degraded  by  a  false  amplification  of  the  noise  floor.  Each  transmit  record  would 
need  to  be  sampled  along  with  measurement  noise.  The  individual  records  would  then  be 
noisier  than  if  they  were  averaged  over  multiple  repetitions  (the  time-domain  equivalent  of  a 
comb  filter).  Here  the  noise  would  represent  the  degree  to  which  the  actual  pulse  shape  and 
timing  could  be  known  rather  than  the  degree  to  which  it  is  consistent  These  two  forms  of 
apparent  displacement  noise  (consistency  of  pulse  vs.  uncertainty  of  pulse)  are  qualitatively 
identical  for  all  practical  purposes 


2.7  Experimental  Data 

The  technique  was  tested  using  an  experimental  setup  that  is  depicted  in  Fig  3.  Shear 
waves  were  excited  in  a  tissue  phantom  ( Blue  Phantom  #  BPF-1410)  by  the  lateral  motion  of  a 
shaker-driven  polycarbonate  sphere  that  was  frictionally  coupled  to  its  surface.  The  shear- 
wave  field  was  then  investigated  using  a  dual-ring  spherically  focused  transducer  (Imasonic 
type  5939)  with  a  2.5  MHz  transmitted  center  frequency.  Displacements  in  and  beyond  the 
phantom  were  examined  to  a  depth  of  36  cm  from  the  transducer's  face  in  the  frequency  range 
of  2  Hz  to  1000  Hz  The  field  depth  and  bandwidth  were  dictated  by  the  spectral  characteristics 
that  were  selected  for  the  transmit  signal,  which  in  turn  were  based  on  properties  of  the 
phantom  (attenuation  and  physical  dimensions).  Processing  the  measured  signals  using  Eq 
(24)  for  a  single  acquisition  (i.e.  V2  second  of  data  and  a  single  position  of  the  ultrasound 
transducers)  yielded  the  displacement  time  waveforms  that  are  depicted  in  Fig  4.  Here,  the 
carrier  signal  was  a  chirp  which  spanned  from  2.1  MHz  to  2.9  MHz  and  the  low-frequency 
excitation  was  a  circular  chirp  that  swept  from  20  to  250  Hz.  The  low-frequency  bandwidth  was 
limited  by  the  shear-wave  source.  The  ultrasonic  transducers  were  2  5  cm  from  the  face  of  the 
phantom,  which  was  submerged  in  water.  This  standoff  was  chosen  to  moderate  the  effects  of 
the  direct  crosstalk  between  the  transmitting  and  receiving  rings.  Weak  scattering  from  the 
water  between  the  transducers  and  the  phantom  is  the  reason  for  the  large  apparent 
displacement  noise  in  this  region,  which  is  depicted  in  frame  A.  Here,  the  noise  can  be  either 
large  or  small  depending  on  the  dominant  source  of  energy  in  the  compressed  carrier  signals 
Where  these  are  dominated  by  the  direct  crosstalk,  the  apparent  noise  is  small  as  the 
denominator  of  the  second  term  in  Eq.  (25)  is  large.  Where  the  compressed  carrier  signal 
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represents  bulk  scattering  from  the  water  or  (more  likely)  small  sidelobes  of  scattering  from 
other  regions,  the  apparent  noise  is  quite  large,  as  the  denominator  in  Eq.  (25)  is  small  In 
neither  case  does  the  noise  represent  actual  displacement.  Inside  the  phantom  (frames  B-H), 
the  apparent  displacement  noise  generally  increases  with  range  because  of  the  attenuation  of 
the  carrier  signal  along  its  two-way  propagation  path.  The  displacement  itself  varies  slowly  and 
shows  dependence  on  range  The  slow  variation  is  because  of  the  low  excitation  frequencies 
involved  At  100  Hz,  the  shear  wavelength  is  about  5  cm  and  the  beam  is  sampling  in  a 
direction  that  is  nominally  across  the  shear  wave  fronts  rather  than  in  the  direction  of  their 
propagation.  The  changes  in  waveform  shape  and  displacement  amplitude  are  because  the 
shear  source  is  directly  above  a  point  midway  along  the  beam,  because  the  angle  between  the 
shear  displacement  and  the  carrier’s  propagation  path  direction  changes  with  range  along  the 
beam,  and  because  the  shear-wave  field  is  composed  of  both  a  direct  signal  and  reflections 
from  the  surfaces  of  the  phantom 

Fig  5  shows  time  delays  extracted  by  correlating  each  measured  displacement  signal  from 
five  different  ultrasonic  transducer  locations  (separated  in  the  X  dimension  in  Fig.  4)  with  the 
drive  signal  for  the  shear  wave  source.  These  delays  show  clear  propagating-wave  structure 
inside  the  phantom,  but  are  primarily  noise-like  in  the  coupling  fluid  outside.  As  would  be 
expected,  the  delay  is  minimal  directly  below  the  shear  wave  source  and  increases  away  from 
it  both  toward  and  away  from  the  ultrasound  transducers.  The  separation  between  the  delays 
for  each  of  the  five  X  locations  is  consistent  with  the  shear  wave  speed  in  the  phantom 
material,  which  was  separately  measured  on  a  section  of  material  by  a  resonant-bar 
technique  XXI  Inside  the  phantom,  the  groups  of  points  that  lay  well  off  of  the  smooth  curve  are 
attributable  to  reflected  shear  waves  from  the  bottom  surface  rather  than  to  noise.  Smaller,  but 
easily  discernable,  correlation  peaks  in  these  records  correspond  to  the  direct  signal.  In  the 
water,  computed  delays  are  not  entirely  noise-like.  Several  of  them  are  consistent  with  an 
extension  of  delays  computed  at  the  back  of  measurement  region.  The  reason  for  this  is  the 
circularity  of  the  carrier  signals.  Since  there  is  very  little  energy  returned  from  the  water,  the 
principle  source  of  modulation  in  this  region  is  wrapped  around  from  later  arrivals  associated 

c  T 

with  scattering  from  the  vicinity  ofz  >  -  —  Here,  the  carrier  may  be  dominated  by  sidelobes  of 
the  pulse-compressed  cross-talk  of  the  transmitted  signal  so  the  amplitude  of  computed 
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displacements  would  not  have  physical  meaning,  but  the  correlation  delays  depend  only  on  the 
phase,  which  is  contributed  entirely  by  the  modulation  and  noise. 

Fig.  6  shows  delays  time  delays  extracted  by  correlating  each  measured  displacement 
signal  from  three  different  ultrasonic  transducer  locations  that  are  separated  in  the  Z  dimension 
in  Fig  4.  Here  it  can  be  seen  that  the  movement  of  the  ultrasonic  transducers  alters  the 
measurements  only  to  the  extent  of  the  corresponding,  and  entirely  expected,  horizontal  shift 
indicating  that  the  post  processing  of  the  data  does  not  introduce  an  unintended  relationship 
between  computed  range  and  time.  An  interesting  feature  of  these  delays  is  the  band 
observed  around  6cm  from  the  first  transducer  location.  This  aligns  with  corresponding 
features  from  the  other  locations  when  the  data  is  shifted  by  twice  the  transducer 
displacement,  as  seen  in  Fig.  4c,  whereas  the  remainder  of  the  data  aligns  when  it  is  shifted  by 
the  transducer  displacement,  as  seen  in  Fig.  4b.  This  feature  is  produced  by  a  second  bounce 
of  the  ultrasound  signal  off  of  the  tissue  phantom  and  the  transducer  face  that  is  larger  than 
the  coincident  speckle  scatter  from  the  interior  of  the  phantom  Thus,  the  translation  between 
arrival  time  and  range  incorrectly  represents  this  portion  of  the  data. 

2.8  Limitations 

Spectral  sampling  limits:  The  temporal  sampling  limit  of  the  displacement  is  defined  by  the 

Nyquist  frequency.  In  the  experiment  this  was  1000  Hz.  There  are  ways  in  which  higher 

frequencies  can  be  sampled  If  the  displacement  has  been  created  by  the  interrogator,  as 
would  be  the  case  in  an  elastography  system,  then  knowledge  and  control  of  its  characteristics 
can  be  used  to  recover  band-limited  signals  above  the  Nyquist  frequency  as  with  any 
conventional  analog-to-digital  converter.xx"  Alternatively,  the  spacing  of  carrier-signal  tones 
may  be  increased  and  the  signal  sampled  at  a  higher  frequency  over  a  shorter  range  if  the 
ambiguity  associated  with  the  resulting  wrap-around  of  signal  could  be  neglected  or  resolved. 
It  is  also  possible  to  select  tones  with  a  larger  separation  in  frequency,  A co ,  with  a  center 
frequency  such  that  the  signal  period  is  not  reduced.  Here  the  cost  would  be  in  spatial 
resolution.  If  the  total  signal  record  is  sufficiently  short,  it  should  also  be  possible  to  process 
the  entire  receive  signal  in  the  frequency  domain  In  this  case,  a  different  spectral  interpolation 
of  the  transmit  and  compression  signals  may  be  used  as  the  modulation  band  associated  with 
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a  particular  tone  will  not  adjoin  that  tone  in  the  frequency,  but  could  still  be  correctly  associated 
in  the  post-processing.  These  sampling  issues  are  similar  in  nature  to  the  problem  of  velocity 
ambiguity  encountered  with  pulse-Doppler  radar  systems. Xl" 

Dynamic  range  limits:  The  requisite  assumption  that  knSm(i) « 1  generally  restricts  the  high 

end  of  the  dynamic  range  for  the  new  technique.  Errors  are  introduced  by  the  inaccuracy  of 
small-argument  approximations  for  functions.  These  errors  limit  the  accuracy  of  the  method  in 
depicting  large  displacements.  The  magnitudes  of  errors  introduced  by  these  approximations 
are  shown  in  Fig.  7.  Based  on  a  0  1  dB  error  tolerance  (40  dB  displacement  signal-to-noise- 
and-distortion  (SiNaD)  ratio),  the  maximum  measurable  displacement  would  be  approximately 
20  mm  with  an  fc=  2.5  MHz  carrier 

Hardware  limitations:  There  are  several  limitations  on  the  technique  that  are  imposed  by  the 
hardware  selected  for  its  implementation.  Chief  among  these  is  the  resolution  of  the  digitizer. 
The  prototype  system  incorporated  a  National  Instruments  N 1-51 05  high-speed  digitizer. 
Although  this  has  a  12-bit  A  to  D  converter,  it  only  delivers  10.5  effective  bits  because  of 
intrinsic  noise  and  distortion. XXI"XXIV  This  introduces  an  artificial  noise  floor  into  all  the  data 
measured  with  the  card.  In  laboratory  measurements  on  a  brain-tissue-mimicking  phantom 
(1.1  dB/cm  attenuation  at  1  MHz),  this  artificial  noise  dominated  the  measured  noise  floor.  The 
N 1-5 1 05  card  also  limits  the  maximum  carrier  frequency.  This  HF  limit  is  primarily  through  the 
speed  with  which  the  card  can  transfer  data  rather  than  its  maximum  sampling  speed  when 
data  is  stored  in  on-board  memory,  as  it  is  capable  of  60  MHz  sampling.  A  second  hardware 
limit  is  imposed  by  the  apparent  spectral  purity  of  the  ultrasonic  transmit  signal.  In  order  to 
properly  implement  the  technique,  the  transmit  signal  must  be  generated  and  recorded  in  such 
a  way  that  it  has  energy  only  at  specific  frequencies.  The  transmitted  energy  in  sidebands 
between  these  tones  must  be  down  by  more  than  100  dB  relative  to  peak  levels  within  a  few 
Hz  of  each  intended  tone  (depending  on  the  lowest  vibration  frequency  that  is  of  interest)  This 
spectral  purity  was  accomplished  by  sharing  a  reference  clock  and  sample  clock  from  an  SRS 
CG-635  clock  source  with  both  the  signal  generator  and  the  digitizer  card.  This  avoided  both 
smearing  the  actual  peaks  of  the  transmitted  spectrum  because  of  clock  jitter  and  smearing  the 
spectral  peaks  of  the  digitized  receive  signal  because  of  clock  mismatch.  The  transmitted 
signals  were  constructed  in  several  different  ways  and  generated  from  a  Tabor  2572A  arbitrary 
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waveform  generator  on  which  records  of  a  single  cycle  of  each  transmit  signal  were  stored  in 
memory  and  continuously  generated.  Data  were  collected  using  periodic  pulses,  linearly  swept 
chirps  and  numerically  optimized  multisinesxxv.  Figure  8  shows  an  example  of  an  experimental 
drive  signal  (a  chirp)  and  its  spectrum  as  it  was  recorded  on  the  N 1-51 05  digitizer.  The  high- 
frequency  roll-off  is  caused  by  an  antialiasing  filter,  as  the  generated  signal  was  flat  across  the 
band  of  interest. 

The  system’s  sensitivity  and  spatial  resolution  are  limited  by  the  characteristics  of  the 
ultrasonic  propagation  in  the  medium  of  interest.  In  order  to  avoid  ambiguity  regarding  the 
position(s)  associated  with  each  displacement  measurement,  it  is  necessary  for  each  full  cycle 
of  the  carrier  to  be  comparable  to  or  greater  than  the  reverberation  time  of  the  sample  under 
study  This  time  is,  in  part,  dictated  by  the  carrier  signal’s  center  frequency  because  of  the 
frequency  dependence  of  attenuation.  Here,  the  limitation  is  that  any  response  to  an  impulsive 
signal  that  occurs  at  t  >T  must  be  significantly  smaller  than  the  response  at  T  >  t  nT  >  0  for 
integer  n  and  that  its  modulation  component  be  similarly  small  This  is  the  same  range 
ambiguity  problem  that  is  encountered  in  pulse-Doppler  radar  systems. Xl"  In  addition  to  the 
amplitude  of  the  low-frequency  vibration,  the  SNR  of  each  target  region’s  measured 
displacement  is  directly  proportional  to  the  energy  backscattered  from  that  region  This  is  a 
function  of  the  scattering  strength  of  the  region,  the  energy  in  the  transmitted  signal,  and  the 
two-way  attenuation  along  the  propagation  path  from  the  source  to  the  region  of  interest.  Thus 
the  presence  of  strong  scatterers  near  the  front  of  a  volume  of  interest  can  enhance  SNR  for 
that  region  while  degrading  it  further  back  in  the  beam.  Attenuation,  while  necessary  to  resolve 
spatial  ambiguity  acts  to  degrade  the  SNR  of  the  measured  displacement  at  longer  ranges  by 
effectively  enhancing  the  noise.  Variability  in  the  strength  of  the  backscatter  throughout  the 
beam  also  occurs  because  of  small-scale  inhomogeneities.  The  manifestations  of  these  effects 
on  the  displacements  computed  using  Eq  (24)  are  that  the  apparent  displacement  noise  varies 
strongly  as  a  function  of  range  from  the  ultrasonic  source  and  receiver. 


2.9  Conclusions 

A  new  technique  for  ultrasonic  vibrometry  has  been  developed  and  tested  The  technique 
offers  several  demonstrable  advantages  over  previously  reported  techniques  for  rapid  full- 
wavefield  measurements  of  small  amplitude  vibrations  In  particular,  the  technique  offers 
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wavefield  measurements  of  small  amplitude  vibrations.  In  particular,  the  technique  offers 
improvements  in  spatial  resolution,  noise  immunity,  and  computational  overhead  to  varying 
degrees  in  comparison  with  either  pure-tone  or  pulsed  vibrometry.  The  constraints  and 
limitations  imposed  by  hardware  on  systems  implementing  the  new  technique  are  generally 
analogous  to  those  encountered  by  similar  systems  implementing  previously  documented 
techniques  or  to  similar  electromagnetic  systems.  In  some  cases,  such  as  the  measurement  of 
relatively  large  displacements,  the  new  technique  is  deficient  in  comparison  to  correlation- 
based  alternatives.  Here,  a  modification  of  the  technique  may  be  used  to  replicate  the  results 
obtained  by  pulsed  vibrometry  without  the  need  to  collect  any  additional  data.  In  doing  this,  it 
should  possible  to  realize  advantages  over  pulsed  systems  by  generating  a  pulse-compressed 
time  record  with  less  noise  than  could  be  directly  achieved  with  pulsed  drive  signals  and  a 
nominally  identical  content  of  underlying  information. 

Development  of  the  technique  is  ongoing  with  current  efforts  focused  on  improvements  to 
system  hardware  and  integration  of  the  vibration  measurement  system  with  an  ultrasonic 
radiation  force  shear-wave  source  for  remote  palpation. 
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List  of  Figures: 

Figure  1:  Periodic  pulse  drive  signal  (A)  one  full  cycle,  (B)  detail  of  pulse  (gray  band 
corresponds  to  a  20-point  post-processing  average),  and  (C)  spectrum  of  a  single  cycle. 

Figure  2:  Comparison  of  apparent  displacement  noise  for  two  different  processors, 
single  pulse  compression  (Eq.  10)  and  the  power-weighted  sum  of  I  and  Q 
compressions  (Eq.  24).  Shown  for  (A)  synthetic  data  and  (B)  measured  data  using  a  1- 
dB  bin  size. 

Figure  3:  Configuration  of  laboratory  ultrasonic  vibrometry  experiment. 

Figure  4:  Measured  displacement  time  waveforms  in  a  tissue  mimicking  phantom 
excited  by  a  surface  mounted  shear  wave  source. 

Figure  5:  Delays  associated  with  correlation  maxima  of  measured  shear-wave 
displacement  with  respect  to  shear  source  drive  signal. 

Figure  6:  Computed  delays  for  five  transducer  locations  relative  to  the  front  face  of  a 
tissue  phantom  (Zo).  (A)  processed  data,  (B)  data  shifted  to  account  for  transducer 
movement,  and  (C)  processed  data  shifted  by  twice  the  transducer  movement. 

Figure  7:  Errors  in  the  computation  of  displacement  associated  with  the  small  argument 
approximations  of  trigonometric  functions  for  f=  2.5  MHz  and  c=1500  m/s. 

Figure  8:  Chirp  transmit  signal:  (A)  full  spectrum  in  the  frequency  domain,  (B)  spectrum 
in  the  near  vicinity  of  a  single  tone  and  (C)  two  cycles  in  the  time  domain. 
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Figure  1 :  Periodic  pulse  drive  signal  (A)  one  full  cycle,  (B)  detail  of  pulse 
(gray  band  indicates  width  corresponding  to  a  20-point  post  processing  average) 
(C)  spectrum  of  a  single  cycle 
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Figure  2:  Comparison  of  noise  for  two  different  processors,  single  pulse  compression 
(Eq,  21 )  and  power-weighted  sum  of  /  &  Q  compressions  (Eq.  31 ),  with 
(A)  synthetic  data  and  (B)  measured  data  ,  plotted  using  a  IdB  bin  size, 
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Figure  3:  Configuration  of  laboratory  ultrasonic  vibrometry  experiment. 


Figure  4:  Measured  displacement  time  waveforms  in  tissue  mimicking  phantom 
excited  by  a  surface  mounted  shear  wave  source 
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Figure  5:  Delays  associated  with  correlation  peaks  of  measured  shear-wave  displacement  with 

respect  to  shear  source  drive  signal 


Figure  6:  Computed  delays  for  five  different  transducer  locations  relative  to  the 
front  face  of  a  tissue  phantom  (ZJ.  (A)  processed  data,  (B)  processed  data  shifted 
to  account  for  transducer  movement,  and  (C)  processed  data  shifted  by  twice  the 
transducer  movement 
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Figure  7:  Errors  in  the  computation  of  displacement  associated  with  the  small- 
argument  approximation  of  trigonometric  functions  for  £=2.5  Mhz,  c=1500  m/s 
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Figure  8:  periodic-chirp  transmit  signal  (A)  full  spectrum  in  the  frequency 
domain, (B)  spectrum  in  the  near  vicinity  of  a  single  tone(2  Hz  bins)  and 
(C)  two  cycles  in  the  time  domain 
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